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' Abstract. The Skyrme-Faddeev-Niemi (SFN) model which is an 0(3) a model in three 

, dimensional space up to fourth-order in the first derivative is regarded as a low-energy 

' effective theory of SU{2) Yang-Mills theory. One can show from the Wilsonian renormali- 

^ zation group argument that the effective action of Yang-Mills theory recovers the SFN in 



the infrared region. However, the theory contains an additional fourth-order term which 
destabilizes the soliton solution. We apply the perturbative treatment to the second deriva- 
tive term in order to exclude (or reduce) the ill behavior of the original action and show 
^ , that the SFN model with the second derivative term possesses soliton solutions. 
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1 Introduction 

The Skyrme-Faddeev-Niemi (SFN) model which is an 0(3) a model in three dimensional space 
^ ' up to fourth-order in the first derivative has topological soliton solutions with torus or knot-like 

structure. The model was initiated in 70's and the interest to it has been growing considerably. 
The numerical simulations were performed in ^121 El IH E] j the integrability was shown in [5], 
^ ■ and the application to the condensed matter physics [Zj and the Weinberg-Salam model |H1 were 

also considered. The recent research especially focuses on the consistency between the SFN and 
fundamental theories such as QCD [9llin | [T] 1 ll2j. In those references, it is claimed that the SFN 
action should be induced from the SU{2) Yang-Mills (YM) action at low energies. One can also 
show from the Wilsonian renormalization group argument that the effective action of Yang- 
Mills theory recovers the SFN in the infrared region However, the derivative expansion for 
slowly varying fields n up to quartic order produces an additional fourth-order term in the SFN 
model, resulting in instability of the soliton solution. 

Similar situations can be seen also in various topological soliton models. In the Skyrme 
model, the chirally invariant Lagrangian with quarks produces fourth order terms after the 
derivative expansion and they destabilize the soliton solution |14l I15j . To recover the stability 
of the skyrmion, the author of JH] introduced a large number of higher order terms in the 
first derivative whose coefficients were determined from those of the Skyrme model by using 
the recursion relations. Alternatively, in Gies pointed out the possibility that the second 
derivative order term can work as a stabilizer for the soliton. 

In this paper, we examine the Gies's supposition by numerical analysis. In Section|21 we give 
an introduction to the Skyrme-Faddeev-Niemi model with its topological property. In Section|31 
we show how to derive the SFN model action from the SU{2) Yang-Mills theory. In Section 
soliton solutions of this truncated YM action are studied. In order to find stable soliton solutions, 
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we introduce a second derivative term which can be derived in a perturbative manner. The naive 
extremization scheme, however, produces the fourth order differential equation and the model 
has no stable soliton solution. In Section |SJ the higher derivative theory and Ostrogradski's 
formulation are reviewed. We show the absence of bound state in the second derivative theory 
using an example in quantum mechanics and introduce the perturbative treatment for the second 
derivative theory. In Sectional we present our numerical results. Section [3 contains concluding 
remarks. 



2 Skyrme— Faddeev— Niemi model 

The Faddeev-Niemi conjecture for the low-energy model of SU{2) Yang-Mills theory is expressed 
by the effective action: 



5sFN = A' / d^x 



where n{x) is a three component vector field normalized as n-n = 1. The mass scale parameter 
A can be scaled out by the replacement Ax — > x and A'^gi — > gi, and for the static energy 
Estt/A Estf Stable soliton solutions exist when gi > 0. 

The static field n{x) maps n : i?^ i— > S*^ and the configurations are classified by the topo- 
logical maps characterized by a topological invariant H called Hopf charge 

H = ^lAAF,F = dA, (2) 

where F is the field strength and can be written as F = {n ■ dn A dn) . 

The static energy from the action Q has a topological lower bound PT] . 

E,tt > KH'/\ (3) 

where K = 4^/2 S^/^vr^^. Note that Ward improved this topological bound by using the Hopf 
map ^Hl- It seems, however, to be an upper bound of the model rather than lower bound. 

Performing numerical simulation, one can find that the static configurations for H = 1,2 
have axial symmetry Thus "the toroidal ansatz" which was studied in [21 is suitable to be 
imposed on these configurations. The ansatz is given by 



ni = y/l - w'^{ri,P) cos{Na + v{t],P)), 
712 = \/l - if^^ (r/, /3) sin(iVa + v{ri, /?)), 

ns = w{r],f3), (4) 

where (r/, /3, a) are toroidal coordinates which are related to the as follows: 

a sinh r] cos a a sinh rj sin a a sin (3 
x = , y = , z = 

T T T 

with T = cosh r] — cos f3. 

The function w^rj, 13) is subject to the boundary conditions w{Q, (3) = 1, ^(oo, [3) = —1 and is 
periodic in (3. v{r], j3) is set to be u(r/, [3) = Mf3 + vo{r], (3) and fo(, j3) is considered as a constant 
map. Equation ((JJ then gives H = NM. 

In this paper we adopt a simpler ansatz than which is defined by 



m = v^l -w'^{r]) cos{Na + M(3), 
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na = ^Jl-w'^{v|) sm.{Na + M(5), 

ng = w{ri), (5) 

where w{rj) satisfies the boundary conditions u;(0) = 1, u;(oo) = —1. We numerically study 
soliton solutions for both ansatz @ and ©. By comparing those results, we find that this 
simple ansatz produces at most 10 % errors and does not much effect to the property of the 
soliton solution. 

By using ©, the static energy is written in terms of the function w{ri) as 

E'stt = 27r^o j dt] Y^^^ + (1 - w'^)UM,N{'n) + ^ sinh77cosh?7(it;')^^M,Af(??) 

The Euler-Lagrange equation of motion is then derived as 

w" ww'"^ , , gi , 2 2 / / \ I 

2+71 +Um,n{ji)w + —^[-2N coth riwUM,N{rj)w 

1 — w-^ (1 — w^y la"^ 

+ ( cosh^ r/ + sinh^ r/) + smhrj coshrjU m,n{v)''J^'') = 0- (6) 

The variation with respect to a produces the equation for variable a. Soliton solutions are 
obtained by solving the equations for a as well as for w. 

3 Effective action in the Yang— Mills theory 
with CFNS decomposition 

In this section, we briefly review how to derive the SEN effective action from the action of SU(2) 
Yang-Mills theory in the infrared limit I13j . Eor the gauge fields A^, the Cho-Faddeev- 
Niemi-Shabanov decomposition is applied |H1 IIUI 1111 112j 

= nC^ + {d^n) xn + W^. (7) 

The first two terms are the "electric" and "magnetic" Abelian connection, and W ^ are chosen 
to be orthogonal to n, i.e. W • n = 0. Obviously, the degrees of freedom on the left- and 
right-hand side of equation ((7J do not match. While the LHS describes 3coior x ^Lorcntz = 12, 

the RHS is comprised of (C^ :)4Lorentz + {n :)2color + {W ^ :)3color X 4Lorcntz - '^n■W^,=Q = 14 

degrees freedom. Shabanov introduced in his paper the following constraint 

x(n,C^,W/,) = 0, with x-n = 0. 
The generating functional of YM theory can be written by using equation as 

Afp and Sg{ are the Eaddeev-Popov determinant and the gauge fixing action respectively, and 
Shabanov introduced another determinant As corresponding to the condition X = 0- YM and 
the gauge fixing action is given by 



'S'ym + Sgi = J d^x 
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Inserting equation ^ into the action, one obtains the form: 



X e 



- (l/2g2 ) J (C^ M^^ C.+Wt, +2W ^, ■ ) 



with 



tW _ a2i 



M, 



K^i^ = di,{n- di,n x (?^n) + d^n ■ d n x n, 
d^{n X 9^n), (in gauge ag = 1) 



w 



and 



Here, has same form of Q^^ but differentiates to the right objects of Q\u- "^^^ 

classical action of n including the gauge fixing term is given by 

5ci = j (fx -^{d^n X d^nf + ^^(^^« x 
The 5 functional is expressed by its Fourier transform 



Integrating over C, W, cf), we finally obtain 

^ e-^^' AfpAs ( det M^) ( det M^) ( det -Qf{M^)~lQf) 



(8) 



where Qf = i{—d^ + d^n ® n) and Qfi differentiates to the right. We perform the derivative 
expansion for the four determinants in equation (jH)) under the following assumptions 

(i) the theory is valid for the momenta p with /c < p < A (/c, A are infrared and ultraviolet 
cut-off), 

(ii) \dn\ < k, 

(iii) the higher derivative terms, such as d'^n are omitted. 



The effective action is then given by 



'S'eff 



(a^n)2 + ^(5^nx5,n)^ + ^^(9, 



2 , 52, 



(9) 



For (7i > and g2 = 0, the action is identical to the FSN effective action (|H). 

In order to get the stable soliton solutions, 52 rnust be positive [^j. However, §2 is found to 
be negative according to the above analysis (see Therefore we consider higher-derivative 

terms and investigate if the model with the higher-derivatives possess soliton solutions. 
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4 Search for the stable soUton solutions (1) 

The static energy is derived from equation @ as 



z o o 



■.= E2{n)+E^l\n) + Ef\n). 



A spatial scaling behavior of the static energy, so called Derrick's scaling argument, can be 
applied to examine the stability of the soliton ^Hl- Considering the map x ^ x' = fix (^u > 0), 
with n(^) = n{fj,x), the static energy scales as 



E2{n) + fi{E^^\n) + E'^>{n)). 



(2), 



(10) 



Derrick's theorem states that if the function e{fi) has no stationary point, the theory has no static 
solutions of the field equation with finite density, other than the vacuum. Conversely, if e(/z) 
has stationary point, the possibility of having finite energy soliton solutions is not excluded. 



equation (fT(i|) is stationary at ^ = y E2/{^E\ + E\ ). Then, the following inequality 

gi{din X djuf + g2{dinf {djuf 

= gi{dinf{djnf - gi{din • djuf + g2{dinf{djnf 
^ g2{din ■ djuf ( •.• {d,nf{djnf ^ (diU ■ djuf) 

ensures the possibility of existence of the stable soliton solutions for g2 ^ 0. As mentioned in 
the Section 13 52 should be negative at least within our derivative expansion analysis of YM 
theory. 

A promising idea to tackle the problem was suggested by Gies ^3]- He considered the 
following type of effective action, accompanying second derivative term 



(11) 



Here we choose positive value of (72 and assign the explicit negative sign to the third term. 
In principle, it is possible to estimate the second derivative term by the derivative expansion 
without neglecting throughout the calculation. 

The static energy of equation with the ansatz (0) is written as 



E. 



stt 



27r a drj 



{w') 



l\2 



92 

4a2 



1 — w"^ 
sinh r] cosh r/ 



+ (1 - w'^)UM,N{'n) + ^ sinh 77 cosh r/(it;')^t/M,Af(r/) 



/\2 



1 — w 



2 + {I - w'^)Um,n{v) 



(w') 

+ (coth ri + sinh^ ri — sinh rj cosh ri) ^ 

1 — 

+ ( sinh cosh - sinh^ rj) (l - w'^)M'^ + 2 < ^ + ^ + ww'Um,n{v) f 

\ [ (1 — ID^)^ 1 — J 



+ sinh r/ cosh ?7 



ww" + (1 - w^)Um,n{v) 
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where w" = The Euler-La grange equation of motion is derived by 



dr/^ \ dw" J dr] \ dw' J dw 
which is too comphcated to write down exphcitly and hence we adopt the foUowing notation 
fo{w,w',w") + gifi{w,w',w")+g2f2{w,w',w",w^^\w'^^^) =0. (12) 

Here w^^\ w^'^^ represent the third and the fourth derivative with respect to i]. The first two 
terms of equation ()12() are identical to those in equation ©. Unfortunately, we could not find 
out stable soliton solutions from equation l\12\\ for any (even in the quite small) value of g2- 
From the relation 

j d'^xKd^n ■ d^n) - (a^n)4] = j d^x{d^n x n)^ 
one easily finds that the static energy obtained from the last two terms in equation (|11|) 

) = f d^x{d'^n X nf (13) 



gives the positive contribution. The total static energy is stationary at fi = y E2/ [E^^^ + E^^'^ 
and hence the possibility of existence of soliton solutions is not excluded. And also, the positivity 
of equation H13() does not spoil the lower bound Q of original SFN and the possibility still 
remains. 

The standard action has a fourth differential order term, but only quadratic in time 
derivative. On the other hand, the action (lllj) involves a term with fourth order in time deriva- 
tive. Then, even if the static energy is constructed only by the positive terms, it is in general 
not bounded from below. As a result, the stability of the soliton is unclear. In the next section, 
we report the basic feature of the higher derivative action and introduce the method to avoid 
the problems inherent by applying an example in quantum mechanics. 



5 Higher derivative theory 

In this section, we make a small detour, i.e., we review the problems in the higher derivative 
theory j2Ul 1211 1221 1231 124j which essentially falls into two categories. The first problem concerns 
the increase in the number of degrees of freedom. For example, if the theory contains second 
derivative terms, the equation of motion becomes up to the order in the fourth derivative. Thus, 
four parameters are required for the initial conditions. If one considers higher-order terms, the 
situation gets worse. However, this is not a serious problem for our study because our concern 
is existence of static soliton solutions. The second problem is that the actions of the theory are 
not bounded from below. This feature makes the higher derivative theories unstable. 

We briefly review the Lagrangian and the Hamiltonian formalism with higher derivative called 
the Ostrogradski method. We consider the Lagrangian containing up to nth order derivatives 

S = j dtC{q,q,...,q^'''^). 

Taking the variation of the action 6S = leads the Euler-Lagrange equation of motion 

t(-'4(|5y)-». 

1=0 
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The Hamiltonian is obtained by introducing n generalized momenta 

j=i+i 

or 

Pn = ^, P^=5^-^^'m, ^ = l,...,n-l, (14) 

and n independent variables 

qi = q, qi = q^''~^\ i = 2,...,n. 

The Lagrangian now depends on the n coordinates qi and on the first derivative (jn = q^^^ ■ The 
Hamiltonian is defined as 

n n—l 

T~(-{H,Pi) = '^PiQi - = '^PiQi+l +Pnqn " ^■ 
i=l 1=1 

The canonical equations of motion turn out to be 

. _dn . _ an 

Qi ~ t; ) Pi — • 
oPi OQi 

Thus, we replace a theory of one coordinate q system obeying 2n-th differential equation with 
a set of 1-st order canonical equations for 2n phase-space variables [qi,Pi]- 

We consider a simple example including second derivative term |24j . defined as 

r -'- /i ( 2 2\ -2 1 2 2 2 "2 

C = -{1 + e UJ )q - ^UJ q - 2^ 9 , 

where constant e works as a coupling constant of second derivative term. The equation of 
motion is 

(l + eV)^' + a;2g + e2g('') = 0. (15) 
From equation (|14)1 . one gets 
dC 



Thus the Hamiltonian becomes 



T-L = q'nq + q-Kq - C 



^ 2 ^ /i , 2 2\ -2 , 2 2 

We introduce the new canonical variables 

wvl — Vl — 



= ,/n ^2..2 (g - ^'?)' = - ^'')' 



VI - e^u^ " ^' eVT 



2,, ,2 
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and the Hamiltonian will have the following form by using these variables 

The corresponding energy spectra is then given by 

{^^l)^' n,m = 0,1,2,.... (16) 

One can see that in the limit e — > the energy goes to negative infinity rather than approaching 
to the harmonic oscillator energy eigenstates. 

To obtain physically meaningful solutions, we employ the perturbative analysis where the 
solution is expanded in terms of the small coupling constant and the Euler-Lagrange equation of 
motion is replaced with the corresponding perturbative equation. The solutions of the equations 
of motion that are ill behaved in the limit e — > are excluded from the very beginning j22|l23lEl] . 

We assume that the solution of equation (fT3)) can be written as 

oo 
n=0 

Substituting equation (|T7j) into equation (fT5|l and taking time derivatives of these equations, we 
obtain the constraints for higher derivative terms 

equation : % + u qo = 0, 

constraints: Qo = —u}'^qo, '4 q = uj"^ go. (18) 

equation : q2 + J^qo + u?q2 + 9 o = 0, 

^ q2 + u?q2 = 0, (using 
constraints : 0.2 = —^"^02, ^ 2 = ^^0.2- (19) 

ry o .... 

equation : 54 + a; g2 + a; 54 + 9 2 = 0, 

^ ^4 + uP'qA, = 0, (using (jlH)), 
constraints : ^4 = —u'^qi, ^ 4 = uj^^q^^. 



Combining these results, we find the perturbative equation of motion up to O(e^) 



which is the equation for harmonic oscillator. 

This perturbative method can successfully exclude the ill behavior of the second derivative 
theory. Of course our problem is different from this quantum mechanics example. While the 
energy spectra ()16|) is quantum value, the energy of the soliton is static and classical. Neverthe- 
less, the perturbative method may give some hint to tackle our problem. In the next section, we 
will apply the method to the action with the second derivative term pip and explore the stable 
soliton solutions. 
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Figure 1. The function w^q) for gi = 0.4, 
g2 = 0, 0.05, 0.1, 0.15, 0.2 (the rescaling radial 
coordinate a; = ry/(l — ry) is used). 



Figure 2. The energy density for gi = 0.4, 
52 = 0, 0.05, 0.1, 0.15, 0.2. 




0.00 0.05 0.10 0.15 0.20 



Figure 3. The energy as a function of g2 {gi = 0.4). 

6 Search for the stable soUton solutions (2) 
(perturbative expansion method) 

As in the case of quantum mechanics, we assume that 52 is relatively small and can be considered 
as a perturbative coupling constant. Thus, in analogy of ((T7j) . the perturbative solution is written 
by a power series in g2 



(20) 



n=0 



Substituting equation H2U|) into equation ()12j) . we obtain the classical field equation in 0{g2) 

fo{wo,WQ,WQ) + gifi{wo,WQ,WQ) = 0. (21) 

Taking derivatives for both sides in equation (|21|) and solving for Wq, we obtain the following 
form of the constraint equations for higher derivatives 

«;« =F»(u;o,w^o), i = 2,3,4. (22) 
The equation in 0(52) can be written as 



(/o + 52/1)0(9!) + f2{wo,w'f),WQ,W^^\w^^^) = 0. 



(23) 
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Substituting the constraint equations (|22() into equation (|23() and eliminating the higher deriva- 
tive terms, one can obtain the perturbative equation of motion 



Now equation ((23]) has stable soliton solutions. 

Our results of the estimated function w{r]) and the energy density are displayed in Figs. ^ 
and [21 (In all figures, we show the results for the case of Hopf charge H = 2; N = 2, M = 1.) 
The dependence of the total energy on g2 is shown in Fig. |21 The change is moderate with 
respect to g2- 

7 Summary 

In this paper we have studied the Skyrme-Faddeev-Niemi model and its extensions by intro- 
ducing the reduction scheme of the SU{2) Yang-Mills theory to the corresponding low-energy 
effective model. The requirement of consistency between the low-energy effective action of the 
YM and the SFN type model leads us to take into account second derivative terms in the action. 
However, we found that such an action including the second derivative terms does not have 
stable soliton solutions. This is due to the absence of the energy bound in higher derivative 
theory. This fact inspired us to employ the perturbative analysis to our effective action. Within 
the perturbative analysis, we were able to obtain the stable soliton solutions. 

Our analysis is based on perturbation and the coupling constant §2 is assumed to be small. 
However, Wilsonian renormalization analysis of YM theory suggests that the coupling con- 
stants gi, g2 (and the mass scale parameter A) depend on the renormalization group time 
t = log/c/A {k,A are infrared, ultraviolet cutoff parameter) and those are almost comparable. 
To improve the analysis, we could perform the next order of perturbation, but it is tedious and 
spoils the simplicity of the SFN model. 

It should be noted that our solutions do not differ much from the solution of original SFN 
model, at least in the perturbative regime. We suspect that an appropriate truncation (such as 
"extra fourth order term -|- second derivative term" ) always supplies the stable solutions that are 
close to the original SFN model. Thus we conclude that the topological soliton model comprised 
of the "kinetic term -|- a special fourth order term" like SFN model is a good approximation. 

Since our results were obtained from numerical study in the perturbative approach, it is 
uncertain whether the existence of the soliton is kept for larger coupling constant g2 ■ To confirm 
that, we should proceed to investigate next order perturbation, or, otherwise, find some analy- 
tical evidence for that. We point out that the perturbative treatment is only used for excluding 
the ill behavior of the second derivative field theory. We hope that applying this prescription 
does not alter the essential feature of the solutions. 

Finally, let us mention the application of the soliton solutions to the glueball. Obviously 
this is one of the main interests to study the model and, many authors have given discussions 
on this subject |13[ I25[ [26] . On the other hand, the possibility of the magnetic condensation 
of the QCD vacuum within the Cho-Faddeev-Niemi-Shabanov decomposed Yang-Mills theory 
have been studied by Kondo [22]. The author claims the existence of nonzero off diagonal gluon 
mass Mx, which is induced in terms of the condensation of the magnetic potential part of the 
decomposition ~ {9fj,n) x n, as 



Throughout our calculation, we scaled out the mass scale parameter A in the action but, in 
this sense, it should reflect the information of such gluon mass, or the condensation property of 
the vacuum. After a careful examination of the value of A, we will be able to make a prediction 
for the glueball mass. 




(24) 



M], = {B^-B,) = {{d^nf). 
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